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Abstract. Virgo is a kilometer- length interferometer for gravitational waves detection located 
near Pisa. Its first science run, VSR1, occured from May to October 2007. The aims of the 
calibration are to measure the detector sensitivity and to reconstruct the time series of the 
gravitational wave strain h(t). 

The absolute length calibration is based on an original non-linear reconstruction of the 
differential arm length variations in free swinging Michelson configurations. It uses the laser 
wavelength as length standard. This method is used to calibrate the frequency dependent 
response of the Virgo mirror actuators and derive the detector in-loop response and sensitivity 
within ~ 5%. 

The principle of the strain reconstruction is highlighted and the h(t) systematic errors are 
estimated. A photon calibrator is used to check the sign of h(t). The reconstructed h(t) during 
VSR1 is valid from 10 Hz up to 10 kHz with systematic errors estimated to 6% in amplitude. 
The phase error is estimated to be 70 mrad below 1.9 kHz and 6 fis above. 



1. Introduction 

The Virgo detector [1] is a kilometer- length interferometer (ITF) located near Pisa (Italy). It is 
designed to search for gravitational waves (GW) in the frequency range from 10 Hz to a few kHz. 
Expected astrophysical sources of detectable GWs are compact objects such as neutrons stars 
or black holes. 



The Virgo first science run (VSR1) was performed from May 18th to October 1st 2007 with 
a sensitivity close to its nominal one, and in coincidence to the end of the fifth science run (S5) 
of the three LIGO detectors [2]. In order to reduce the false detection rate and improve the 
estimation of the parameters of potential sources, the data of all the detectors are used together 
to search for a GW signal. 

The purpose of the Virgo calibration is to measure the parameters needed (i) to estimate 
the ITF sensitivity to GW strain as a function of frequency, h(f), and (ii) to reconstruct from 
the ITF data the amplitude h(t) of the GW strain signal. To achieve optimum sensitivity, the 
positions of the different mirrors are controlled [3] to have, in particular, beam resonance in the 
cavities, destructive interference at the ITF output port and to compensate for environmental 
noise. The controls would also attenuate the effect of a passing GW below a few hundreds hertz. 
A synopsis of the longitudinal control loop and its components is given in figure 1. The effects 
of the longitudinal controls have thus to be precisely calibrated in the frequency domain. Above 
a few hundreds hertz, the mirrors behave as free falling masses in the longitudinal direction. 
The main effect of a passing GW would then be a frequency-dependent variation of the power 
at the ITF output, characterized by the ITF optical response. Finally, the readout electronics 
of the output power and its timing precision have also to be calibrated. 

In this paper, the main method used to get an absolute length calibration is first described 
and the mirror actuator calibration results are highlighted. The way the Virgo sensitivity is 
measured is then shown. In the last section, the principle of the h(t) reconstruction is described 
as well as some measurements performed to estimate the errors on the h(t) signal. 
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Figure 1. Synopsis of the Virgo longitudinal control loop. A variation of the differential length 
of the arms of the ITF, AL, is sensed through the optical response which converts the length 
variation into laser power variation at the ITF output port. The power impinging on photodiodes 
is read and sampled at 20 kHz in the channel Pacp (W). This channel is used as error signal to 
compute the control signals (V) sent to the different mirror actuators. 

Calibration signals can be injected into the loop to force differential length variations AL ca nf)- 
The standard path is through the electromagnetic actuators. An auxiliary setup, the photon 
calibrator (peal), can also be used (see section 4.3). 



2. Absolute length measurement for mirror actuation calibration 

The ITF calibration is based on absolute length measurements. The displacement induced by 
the mirror actuators is reconstructed using the ITF as a simple Michelson. The non-linear 
method used to reconstruct the arm differential motion in the freely swinging Michelson with 
passing fringes is described. The laser wavelength is used as length standard. 

Mirrors of the Virgo ITF are mis-aligned to get a simple Michelson configuration. Different 
configurations are used to calibrate the beam-splitter and arm mirrors actuations (comprising 
the beam-splitter, one end mirror and the input mirror in the other arm). The typical results 
of the mirror actuation calibration are given. 



2.1. Absolute length measurement 

In a simple Michelson ITF, the phase difference A$ between the two interfering beams is 
function of the differential arm length AL, using the laser wavelength A = 1064 nm as standard: 
A*(t) = f AL(t). 

In a simple Michelson ITF with frontal phase-modulation, the power of the output beam (Vdc) 
as well as the demodulated power (Vac) signals are functions of the phase difference A<I> between 
the two interfering beams: Vdc = — 7cos(A<I>)) and Vac = asin(A<3?), where a and (3 
are proportional to the laser power and 7 is proportional to the ITF contrast. Therefore, in the 
(Vdc,Vac) plane, the signal follows an ellipse as shown in figure 2. 

The measurement of the differential arm length AL requires a non-linear reconstruction. 
The ellipse is fitted using the method described in [4]. The fit gives the ellipse center position 
(theoretically (/3,0)) and the axis widths (theoretically (f3j, a)). The variations of the parameters 
a, f3 and 7 are monitored and found to be of the order of 0.5% during a few-minute dataset. 
A^' can then be estimated directly from the ITF signals as the angle between the ellipse axis 
along DC and the line from the ellipse center to the current point position (Vdc^Vac)- Using 
a suitable ellipse tour counting, the right number of 2tt is added to A$' to recover completely 
the angle A<I>. The differential arm length AL rec is then computed from equation 2.1. 
Figure 3 is an illustration of the method. In the window At, 6 interfringes passed on the DC 
signal: It indicates a differential arm elongation of 6 x -| = 3.19 /mi. In the same window, the 
reconstructed AL varies by ~ 3.18/iin. 
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Figure 2. Vac vs Vdc measured in free 
swinging Michelson data (configuration 
comprising the beam-splitter and the two 
arm input mirrors). 
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Figure 3. Red curve (top): measured time 
series of Vdc in free swinging Michelson data, 
in arbitrary units. Black curve (bottom): 
reconstructed time series of AL, in meters. 



2.2. Mirror actuator calibration 

The Virgo electromagnetic actuators are used to induce a longitudinal motion to the suspended 
mirrors from a voltage signal. As shown in figure 4, the signals pass through some electronics 
and pairs of coil- magnet (^4) and the pendulum mechanical response (P). Their calibration 
consists in measuring the transfer function between the input voltage V{ n j and the induced 
motion ALj nj -. This is done setting the ITF in free swinging Michelson configurations: sine 
signals Vinj are injected to the actuators and the motion AL rec is resconstructed using the 
method described in section 2.1. The complete response of the ITF is shown in figure 4: the 
reconstructed AL rec has to be corrected for some effects to get the true induced mirror motion 
ALi n j, mainly some delays due to the light propagation time in the arms (O) and due to the 
output power readout electronics (S). 



Thirteen sine signals have been injected from 5 Hz to 1. kHz. Such injections were done 
at different amplitudes in order to check the linearity of the measurements. They were also 
performed every two weeks during VSR1 in order to monitor their time stability. The mirror 
actuation responses ALj n j jV% n j are measured with a statistical error better than 1%. A regular 
monitoring of the actuation has shown ~ 2% variations and different measurements had up to 
~ 2% offsets from the standard measurements. Systematic errors of the order of 4% on the 
modulus and 50 mrad on the phase have been derived. 
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Figure 4. Signal processing in the free swinging Michelson data used for the mirror actuation 
calibration. See text for details. 



3. Virgo ITF response and sensitivity measurements 

The first goal of the calibration is to measure the Virgo sensitivity. It is based on the knowledge 
of a mirror actuator response calibrated as shown in section 2. The ITF response in closed-loop 
is first measured and then used to estimated the ITF sensitivity. 

3.1. Virgo ITF response 

The ITF response, in W/m, is the the transfer function from the differential arm length variation 
AL to the ITF output power Vac?- When the Virgo ITF is in closed- loop (see figure 1), some 
white noise Vi n j is injected up to 1.5 kHz through an end mirror actuator. Using the actuator 
calibration, the induced ALj n j is derived. The transfer function AL{ n j /Vacp (inverse of the 
ITF response) is then measured up to ~ 1.5 kHz. 

The longitudinal controls do not act above a few hundred Hz: the ITF response is only the 
optical response and the output power sensing. A simple model is thus fit to the measured 
transfer function from 800 Hz to 1 kHz. The optical response is modeled by a simple pole at 
500 Hz due to the arm cavities with finesse 50. The finesse of the arm cavities changes over time 
by as much as 3%. These variations are ignored at this stage. The gain of the optical response 
is let free in the fit. The measurements of the output power sensing are not described in this 
paper: the fitted model of the sensing lies within 2% in modulus and 20 mrad in phase from the 
measured response up to 6 kHz. 

An example of a measured inverse ITF response during VSR1 is given in figure 5. 

3.2. Virgo sensitivity 

The Virgo sensitivity is measured from the ITF response in closed-loop, in m/W, and from the 
spectrum of the output signal Vacp> m W: they are multiplied to get the sensitivity in meters. 
It is then divided by the Virgo arm length, 3 km, to get the sensitivity in strain as shown in 
figure 6. The sensitivity is measured a few minutes after the ITF response in order to avoid 
possible variations of the optical gain. 

Note that the sensitivity is measured with systematic errors of the order of 4% coming from 
the mirror actuator calibration below 1 kHz. The errors above 1 kHz are higher (5-10%) due 
to (i) the use of a fixed cavity finesse, (ii) the power sensing model errors and (iii) the fact that 
the controls were not totally negligible in the fitted frequency range during VSR1. 




Figure 5. Virgo ITF response measured in 
September 2007. Black: measured transfer 
function up to 1.5 kHz (with increasing 
statistical uncertainties above a few 100's Hz). 
Red: extrapolated transfer function. 



Figure 6. Virgo sensitivity measured in 
September 2007. 



4. h(t) strain reconstruction 

The second aim of the calibration process is the reconstruction of the GW strain time series 
h(t). The principle of the reconstruction used for the VSR1 data is described. Some estimations 
of the errors and the sign of h(t) are then given. 

4-1. Reconstruction principle 

The principle of the reconstruction of h(t) from the ITF output power and the mirror longitudinal 
controls is given in figure 7. The effective arm differential length variation AL e ff ect i ve is 
reconstructed in the frequency domain from the output power Vacp and the calibrated optical 
response and power sensing. It contains the contribution from the noise, from the controls 
and possibly from the GWs. The contributions from the controls are reconstructed in the 
frequency domain from the correction signals sent to the mirror actuators knowing the actuator 
responses. They are then subtracted from AL e ff ec ti ve to get the final AL with the noise and 
GW contributions only. The result is then divided by 3 km and put back in the time domain 
to get the strain h(t). The contribution from the power line (50 Hz and harmonics) is finally 
subtracted. 

4-2. Error estimation using mirror actuators 

In order to estimate the errors on the reconstructed strain h(t), some datasets with sine signals 
Vi n j injected to an out-of-loop mirror (an input mirror) were taken every week during VSR1. 
It corresponds to an injected strain hi n j = ALj nj / (3000 m). The transfer function of the 
reconstructed over injected strain h/hi n j is then measured. This transfer function is expected 
to be flat: at in phase and at 0.97 in modulus (since the motion of the input mirror changes 
both the arm length and the Michelson length, the optical response of Virgo is 3% lower for an 
input mirror than for an end mirror). 

The differential arm length variations ALj n j induced by the injected signal is derived from 
the mirror actuation calibration, with errors of the order of 4% in modulus and 50 mrad in 
phase. 

The transfer function measured every week during VSR1 is shown in figure 8. The observed 




Figure 7. Synopsis of the reconstruc- 
tion of h(t) from the ITF output power 
"Pacp and the longitudinal correction sig- 
nals sent to the mirror actuators {zCorr 
channels). See text for details. 



Figure 8. Top curve: ratio between the 
reconstructed amplitude of the lines in hit) and 
the expected injected amplitude. Bottom curve: 
phase difference for the two signals. The difference 
symbols correspond to different measurements 
performed once per week during VSR1. 



time variations at a given frequency are compatible with statistical errors. The difference 
between the measurements and the expected values gives an estimation of the systematic errors 
on the reconstructed strain hit). They are of the order of 4% in modulus and 50 mrad in phase 
below 1 kHz. This is compatible with the estimated systematic errors of the mirror actuation 
calibration used in the h(t) reconstruction. 

4-3. "Photon calibrator" and sign of h(t) 

As in GEO [5] and LIGO [6, 7], a photon calibrator (peal) was setup in Virgo during VSR1. It 
is used as an independent mirror actuator: the radiation pressure of an auxiliary laser is used to 
displace one input mirror. Its main use during VSR1 has been to check the sign of hit), defined 
in common with LIGO as h = x L y where, in Virgo, x and y are the north and west arms 
respectively. 

Sine signals injections were performed with the Virgo peal, corresponding to an expected 
signal hi n j. The phase between the reconstructed h and the expected one hi n j has been measured 
below 1 kHz. It is expected to be 0. The phase at low frequency converges to within 15 mrad: 
this proves that the sign of the reconstructed h(t) is correct. 

Note that the peal has been used also to cross-check the mirror actuation standard calibration. 
It agrees within the large peal systematic errors of 20%. 

5. Conclusion 

The Virgo standard calibration is based on a non-linear reconstruction of the error signal in free 
swinging Michelson configurations. It relies on the laser wavelength as length standard. It has 
been used to calibrate the mirror electromagnetic actuators within 4% in modulus and 50 mrad 
in phase. 

The Virgo photon calibrator has been mainly used to check the sign of h(t). 

The principle of the reconstruction of the strain h(t) applied on the VSR1 data has been 
given along with some estimations of the errors. Some other sources of errors must be taken into 
account, especially from the ITF output power sensing and from some timing issues. Taking 
them into account, the reconstructed h(t) during VSR1 is valid from 10 Hz up to the Nyquist 



frequency of the channel (2048 Hz, 8192 Hz or 10000 Hz) with systematic errors estimated to 
6% in amplitude. The phase error is estimated to be 70 mrad below 1.9 kHz and 6 fis above. 
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